Astronomy & Astrophysics manuscript no. Ibigot 


©ESO 2011 


October 6, 2011 





Letter to the Editor 



The diameter of the CoRoT target HD 49933 
Combining tlie 3D iimb darl<ening, asteroseismology, and interferometry 

L. Bigot^ D. Mourai-d^, P. Berio^, F. Thevenin^ R. Ligi^, I. Tallon-Bosc^ O. Chesneau^, O. Delaa^, N. Nai'detto^, K. 
Perrauf*, Ph. Stee^, T. Boyajian^ P Morel\ B. Pichon\ P Kervella'', F. X. Schmider^, H. McAlister^•^ T. 
ten Brummelaar^, S. T. Ridgway^, J. Sturmann^, L. Sturmann^, N. Turner^, C. Fanington^ and P. J. Goldfinger^ 



o 

(N 

o 
O 
in 



o 



> 

in 

00 

o 



X 



Universite Nice-Sophia Antipolis, Observatoire de la Cote d' Azur, CNRS UMR 6202, BP 4229, F-06304 Nice Cedex, France. 
Universite Nice-Sophia Antipolis, Observatoire de la Cote d' Azur, CNRS UMR 6525, BP 4229, F-06304 Nice Cedex, France. 
UCBL/CNRS CRAL, 9 avenue Charles Andre, 69561 Saint Genis Laval cedex, France 

UJF-Grenoble 1 / CNRS-INSU, Institut de Planetologie et d'Astrophysique de Grenoble, UMR 5274, Grenoble, F-38041, France. 
CHARA and Department of Physics and Astronomy, Georgia State University, P.O. Box 4106, Atlanta, GA 30302-4106, USA 
LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Universite Paris Diderot, 5 place Jules Janssen, 92195, Meudon, France 
Georgia State University, PO. Box 3969, Atlanta GA 30302-3969, USA 
CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA 
National Optical Astronomy Observatory, P. O. Box 26732, Tucson, AZ 85726, USA 



Received ; accepted 



ABSTRACT 



Context. The interpretation of stellar pulsations in terms of internal structure depends on the knowledge of the fundamental stellar 
parameters. Long-base interferometers permit us to determine very accurate stellar radii, which are independent constraints for stellar 
models that help us to locate the star in the HR diagram. 

Aims. Using a direct interferometric determination of the angular diameter and advanced three-dimensional (3D) modeling, we derive 
the radius of the CoRoT target HD49933 and reduce the global stellar parameter space compatible with seismic data. 
Methods. The VEGA/CHARA spectro-interferometer is used to measure the angular diameter of the star A 3D radiative hydrody- 
namical simulation of the surface is performed to compute the limb darkening and derive a reliable diameter from visibility curves. 
The other fundamental stellar parameters (mass, age, and T^if) are found by fitting the large and small p-mode frequency separations 
using a stellar evolution model that includes microscopic diffusion. 

Results. We obtain a limb-darkened angular diameter of ^ld = 0.445 ± 0.012 mas. With the Hipparcos parallax, we obtain a radius 
of R = 1.42 ± 0.04 Rq. The corresponding stellar evolution model that fits both large and small frequency separations has a mass of 
1.20 ± 0.08 Mo and an age of 2.7 Gy The atmospheric parameters are T^g = 6640 ± 100 K, logg= 4.21 ± 0.14, and [Fe/H] = -0.38. 

Key words. Stars: Fundamental parameters. Oscillations - Methods: Numerical - Physical data and processes : Asteroseismology, 
Convection, Hydrodynamics, Radiative transfer - Techniques: Interferometric - Individual: HD49933 



1. Introduction 

The F5V solar-like pulsator HD 49933 (HR 2530, HIP 32851) 
is one of the prim ary targets of the CoRoT space mission 
(iBa^hn et al.ll2006l) . The discovery of its oscillation was made 
fro m the ground by obs erving highly resolved time-series spec- 
tra (iMosser et al.ll2005h . Since the launch of the spacecraft in 
December 2006, two runs of 60 and 137 days have been per- 
formed leading to a rich p-mode spectrum of 5 1 frequencies cor- 
responding to the { — 0, 1,2 eigenmodes ('Appou rchaux et all 
[2008; Benomar et al ] l2009l: fKallingeret al. 2010). 

The interpretation of these eigenmodes to derive con- 
straints on the stellar interior depends on the accuracy of 
the determination of the fundamental parameters of the star 
(Tetf, log [Fe/H]). All chemical abundance analyses of HD 
49933 converge to that of a slightly metal-poor star with JFe/H ] 
between -0.3 and -0.5 (e.g. lGillon & Magainl2006tlBrunttl200"9l) . 
The determination of its effective temperatu re and gravity have 
been revisited several times. As reported in iBrunttf (l2009h . the 



Teff obtained by different methods ranges from 6450 to 6780 K, 
which is by far too uncertain to tightly constrain the stellar evolu- 
tion models in the HR diagram. Similarly, log g is also not very 
well constrained with values ranging from 4.0 to 4.3. Several 
authors have attempted to determine Teff and logg by fitting 
the Hfl and Mg i lin es, respectively dRvabchikova et all 120091: 



lKalIingeretal.ll201flh . Their solutions are most consistent with a 
cooler temperature (^ 6500 K) and a lower gravity (w 4.0). 
Since spectrometry and photometry lead to stellar parameters 
with a large uncertainty, in this work we use a different approach 
by determining its diameter using interferometry. Combined 
with seismology, it p rovides a very good constraint of the stel- 
lar parameter space (ICreeve v et al.' '2007*). This approach was 
proposed for the first time by Kervella et al. (2003) for the bi- 
nary star ffCen A&B and later wi th similar success for other 
asteroseismic targ ets ('Cun ha et al.| 
iBruntt et al.ll2009l: IMazumd ar et aU 
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In this letter we present the first campaign of interferometric 
meas urements of HP 49933 obtained using the VEGA instru- 
ment dMourard et al.ll2009HMourard et al.ll201 ih at the CHARA 
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Fig. 1. Monochromatic (730 nm) center-to-limb emergent intensity at a representative time of the simulation for various cosines 
of the viewing angle : ju = 1.0 (disk center) down to 0.1 (towards the limb). The horizontal sizes of the simulation domain are 
21000 x21000 km. 



Array dten Brummelaar et al.l2005h . We used state-of-art 3D hy- 
drodynamical simulations of the surface of the star to calculate 
the limb darkening and to compute the visibility curves prior ex- 
tracting the angular diameter The derived stellar radius is then 
used to constrain the other fundamental stellar parameters by fit- 
ting the oscillation frequency separations. 

2. CHARA/VEGA observations 

Observations of HD 49933 were performed on October 16, 2010 
using the VEG A instrument of the CH ARA Array and the instru- 
ment CLIMB dSturmann et al.ll2010l) for 3T group delay track- 
ing. The telescopes El, E2 and W2 were used giving access to 
ground baselines of 66, 156, and 221m. The seeing was stable 
with a value of rO between 10 and 13 cm. For the absolute cal- 
ibration of the squared visibilities, we observed two different 
calibrators, CI = HD 55185 and C2 = HD 46487. The observ- 
ing sequence was C1-T-C2-T-C2-T-C1, each block being almost 
20mn. The data were processed using the standard procedur e 
of the VEGA instrument as described in iMourard et al.l (1201 ih . 
To validate the absolute calibration, we first considered CI as a 
target and thus estimated its uniform disk (UD) diameter using 
C2 as the calibrator. A spectral band of AA = 30 nm centered 
around A = 735 nm was used for the data processing. The equiv- 
alent UD diameter of C2 is taken as 0hd46487 = 0. 180+0.013 mas 
accor ding to the SearchCal tooQ of the JMMC dBonneau et alJ 
l2006h . Fitting a UD model to the three squared visibilities pro- 
vides an estimate of the angular diameter of HD 55185 of 
%D55i85 = 0.474 + 0.014 mas. Spectral bandwidth and baseline 
smearing are included but only introduce a bias less than 0.1% 
on the diameter determination. Surface brightness relations, as 
provided by SearchCal, provide an independent estimate of the 
angular diameter We found 0.478 mas and 0.476 mas from the 
(V - R) and (B - V) relations, respectively. Our result agrees 
very well with these photometric estimations. We finally use CI 
and C2 to calibrate the measurements on the target. To avoid a 
strong telluric feature, data are processed using a spectral band 
of AA = 20 nm centered around A = 740 nm. The final error 
bars in the target measurements take into account the measure- 
ment noise, the dispersion in the individual measurements, the 
diameter of the calibrators, and their errors. 

3. Determination of the diameter 

To derive a reliable diameter of HD 49933, we base our analy- 
sis on a 3D hydrodynamical simulation of the surface of the star 
[Allende Prieto et al. ( 2002) demonstrated that for Procyon, a star 
similar to HD 49933, the 3D limb darkening is significantly dif- 
ferent from the corresponding one-dimensional (ID) hydrostatic 

' Available at ,http ://w w w.jmmc .fr/searchcal | 



and homogeneous one. The differences were found to be larger 
in the visible (the VEGA domain) than the infrare d. We adopt 
the sa me methodol ogy as used for both a Cen B dBigot et al.l 
l2006h and Procyon dAufdenberg et al.ll2005h . A hydrodynamical 
simulation of the surface is used to obtain 3D time-dependent 
limb-darkened intensities, instead of the simplified ID homoge- 
neous or uniform disk models, to compute the visibility curves. 
The angular diameter is found by minimizing the differences be- 
tween the synthetic and observed visibility curves. 

3.1. The 3D simulation of the surface 

We use the state-of-the-art radiative hydrodynamical code 
(stagger code, Nordlund & GalsgaarcQ) to simulate the surface 
convection and stratification of HD 49933. In a local box, the 
code solves the full set of conservative hydrodynamical equa- 
tions coupled to an accurate treatment of the radiative trans- 
fer. The code is based on a sixth order explicit finite differ- 
ence scheme. The equations are solved on a staggered mesh 
where the thermodynamical variables are cell centered, while 
the fluxes are shifted to the cell edge. The domain of simula- 
tion contains the entropy minimum located at the surface and is 
extended deep enough to have a flat entropy profile at the bottom 
(adiabatic regime). The code uses periodic boundary conditions 
horizontally and open boundaries vertically. At the bottom of 
the simulation, the inflows have constant entropy and pressure. 
The outflows are not constrained and are free to pass through 
the boundary. We used a realistic equation-of-stat e that accounts 
for io nization, recombination, and dis sociation dMihalas et al.l 
1 19881) and continuous + line opacities dGustafsson et al.l 20081) . 
Radiative transfer is crucial si nce it d rives convection through 
entropy losses at the surface dStein & N ordlund 19981) and is 
solved using the Feautrier's scheme along several inclined rays 
(one vertical, eight inclined) through each grid point. The wave- 
length dependence of the radiative transfer is taken into account 
using a binning scheme in which the monochromatic lines are 
collected into 12 bins. The numerical resolution used for the 
present simulation is 240"^. The geometrical sizes are 21 x 21 
Mm horizontally and 5 Mm vertically. The horizontal sizes of 
the domain are defined to contain a suflicient number of gran- 
ules at each time-step and the vertical one to be deep enough to 
ensure an adiabatic regime at the bottom. The stellar parameters 
fliat define our 3D model are Teft = 6690 ± 35 K, logg = 4.21 
and a scaled solar chemical composition dAsplund et al.ll2009l) 
down to -0.5 dex. The uncertainty in Tetf represents the fluctua- 
tions with time around the mean value. The simulation was run 
for a couple of stellar hours to get the full hydrodynamical re- 
laxation. The convection in such an F star is more efficient than 



^ 1995, |http://www.astro.ku.dk/~kg/Papers/MHDcode.ps.gz| 
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Fig. 2. Our best-fit model of the observed squared visibilities 
(black dots) with the calculated one (full line) with a reduced 
= 0.47. The angular diameter is 0ld - 0.445 mas. 



for the Sun: the rms vertical velocity is 3.9 km/s, which is about 
twice the solar value. The granulation pattern is shown in Fig.[T] 



3.2. Limb darkening and visibility curves 

Hereaf ter, we use the methodology presented in iBigot et al] 
(l2006l) . The snapshots are taken after relaxation and cover a se- 
quence of several convective turnover times. We extract snap- 
shots every five minutes over a sequence of one hour. For each 
snapshot, we calculate the radiative transfer in the selected wave- 
length bandwidth of the VEGA instrument AA = [730 - 750] 
nm to get the monochromatic limb-darkened intensities for each 
time step, each position at the stellar surface, and each cosine 
of the angle made between the vertical and the line of sight. 

The normalized fringe visibilities are obtained using the 
Cittert-Zernike theorem 



< Ia(m) > Jo {ttBOA-' ^Jl-fi^)ndn 
Jo < Ia{m)> l^dju 



(1) 



where B is the projected baseUne (in meters), Jq is the zero order 
of the Bessel function, and < IaQj) > is the time- and horizontal- 
average of the 3D intensities. We integrate the monochro- 
matic visibilities over the spectral domain AA of the instrument 
weighted by the spectral transmission T(A) 



J^^[VAB,e)TiA)fdA 



LATiDfdA 



(2) 



The function T{A) is derived from the direct flux measurement 
of the target and checked on several bright stars. 

3.3. Thex^ fit of the data 

The angular diameter is determined by a Levenberg-Marquardt 
least-squares fit to the observed squared visibilities with the 
model. The result is shown in Fig. |2] The value derived from 
the 3D limb darkening is ^ld = 0.445 + 0.012 mas. The re- 
duced of the fit is 0.47. U sing the revised Hip parcos par- 
allax n - 33.69 + 0.42 mas (Ivan Leeuwenll20()7h . we obtain 
a linear radius of R - 1 .42 + 0.04 Rq. The corresponding ra- 
dius derived from ID hydrostatic limb-darkening (Claret 200^ 
is ^ - 1 .45 + 0.04 Rq, which is significantly larger by 1 cr be- 
cause of the more pronounced limb darkening in ID. We note 



that the derived radius is not very sensitive to the choice of stellar 
parameters (those selected in Sect l3.1l ). We indeed tried several 
values of Teff (= 6500, 6750 K) and log g (= 4.0, 4.5) and found 
very small differences in terms of derived radius (< 0.002 /?©). 
The derived angular diameter 0ld agrees with the prediction of 
0.452 + 0.007 mas ob tained by the surface brightness relation 
(iThevenin et al.ll2"006h . The relative uncertainty in the angular 
diameter (2.7%) is larger by an order of magnitude than the 
smallest uncertainty ever achieved for an interfer ometric radius, 
(0.2% ) derived in the case of the a Centauri A dKervella et al.l 
I200I . However, we emphasize that the target in the present 
work is 20 times smaller This uncertainty dominates the to- 
tal uncertainty of about 3%, including the error in the parallax. 
Nonetheless, the relative uncertainty in the radius is sufficiently 
accurate to tightly constrain the mass using asteroseismic fre- 
quencies dCreevev et al.ll200'7h . 

4. Asteroseismic fits of acoustic frequencies 

We estimate the mass and the age of HD49933 using astero- 
seismic data. The aim is to find a stellar evolutionary model 
whose radius equals the interferometric one , whose Teff, log;?, 
and [Fe/H] match the spectrometric ones (lBrunttll2009h . and 
whose eigenmode fre quency separations match the CoRoT data 
(iBenomar et al.ll2009h . We compute a grid of models of mass be- 
tween 1.7 and 2.1 M©, and search for the one that most closely 
matches the observed small and large frequency separations. 
The stellar model evolves until the radius equals the interfer- 
ometric one. The initial metallicity is adjusted to match the 
spectroscopic [Fe/H]. The stell ar evolutionary models are com- 
puted using the CESAM2k code (lMorellll997t iMorel & LebretonI 



I2008h . In addition, we adopt Chemical diffusion dBurgerslll969h 
combined with extra-mixing dMorel & Theveni nl l2002h and the 
CanujB.&.Mjzzitelli ( 199-L) formalism of convection. The ro- 
tation of HD 49 933 has a minor impact on the modeling 
dPiau et al. I l2009h and is ignored. For each stellar model, we 
compute eigenfrequencies between 11 00 and 2500 fiViz for dif- 
ferent Legendre polynomial degrees I - 0, 1 , 2 by solving 
the standard linearized equations of continuity and m omen- 
tum for adiabatic oscillations (see e.g. lUnno et aDll989l) . Since 
these modes pulsate in high overtones (n), we added the shift 
in frequency caused by the c entrifugal distortion of the star 
dDziembowski & Goodell 19921) . 

We fit both the large and small frequency separations that we 
define as follows 

A„,f = ^ (v„+i.f - v„_i,f) and 6„,( = v„/ - ^ (v„+i.f + v„,c+i) (3) 

Since the radius is given, the large separations are mainly sen- 
sitive to the mass of the star The small separations are more 
sensitive to the deep interior and are therefore a good indicator 
of the stellar age (e.g.;Roxburgh & Vorontsov 2003). 

The solution that minimizes the differences between the ob- 
served and computed separations is given in Tab. [1] and the 
corresponding frequency separations are shown in Fig. [3] We 
note that a core overshoot of 0.35 is needed to fit the small 
separations. We get a better fit of the small separations by in- 
cluding the centrifugal effects on the oscillations. The fits to 
both smal l and large separation s are very satisfactory and bet- 
ter than in lBenomar et al. To roughly estimate the errors 
in the mass and Teff, we use the observed errors in the average 
large separation (~ M^^^R^^I^), in the luminosity L (~ R^T^^), 



Code kindly provided by W. Dziembowski 
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Table 1. Our stellar evolution model for HD 49933. The mass M, initial helium content Yq, metallicity (Z/X)o, core overshoot ffov, 
and mixing length a are adjusted to reproduce the radius obtained by interferometry and the observed large and small separations. 
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Fig. 3. Large (left panel) and small (right p anel) separation s for { - (blue) and ( - I (red) calculated with the stellar model of 
table [T] The observed data (bullets) are from lBenomar et all (l2009l) . For each frequency, we took the maximum error bar at Icr. 



and in the interferometric radius. Ou r stellar evolution mode l 
agrees well with one of the solutions of lCreevev & Bazoi(l201 ih . 
who used a local minimization method based on radial modes. 
The mass (1.20 + 0.08 Mp) is close to that of other works 
(lAppourchaux et al ] l2008l: iBenomar et al.ll20Tl iKalUnger et all 
I201(]h . The luminosity log L/Lp - 0.55 matches the one ob- 
served (iMichel et al.l2008l) . The derived atmospheric parameters 
(Teff = 6640 + 100 K, logg= 4.21 + 0.14, [Fe/H] = -0.38) are 
in good agreement wi th the values obtained by Stromgren pho- 
tometry (lBrunttll2009l). Our T^ff is larger than the value 6500 
K) de rived from (iRvabchikova et al.1 l2009t iKalUnger et"aD 
l2oToh . However, we emphasize that this strong line is not an ex- 
cellent temperature indicator since it spreads over a wide wave- 
length range. 



5. Conclusion 

We have reported our long-baseline interferometric observations 
of the CoRoT target HD 49933 using the VEGA instrument at 
the CHARA Array. Combined with 3D hydrodynamical simula- 
tions of the surface, we derived its radius, R = 1.42 + 0.04 R©, 
which agrees with the radius derived from the surface brightness 
relation. The use of this interferometric radius reduces the stel- 
lar parameter space and helped us to more tightly constrain the 
position of the star in the HR diagram. The fit of the p-mode fre- 
quencies led us to determine the other stellar parameters. The de- 
rived atmospheric ones (Teg - 6640 + 100 K, logg- 4.21+0.14, 
[Fe/H] = -0.38) agree well with photometric studies. 
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